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We  report  on  a  study  of  the  electronic  and  magnetic  properties  of  the  triangular  antiferromagnetic 
{CU3}  single-molecule  magnet,  based  on  spin  density  functional  theory.  Our  calculations  show  that 
the  low-energy  magnetic  properties  are  correctly  described  by  an  effective  three-site  spin  s  =  1/2 
Heisenberg  model,  with  an  antiferromagnetic  exchange  coupling  J  ft  5  meV.  The  ground  state 
manifold  of  the  model  is  composed  of  two  degenerate  spin  S  =  1/2  doublets  of  opposite  chirality. 

Due  to  lack  of  inversion  symmetry  in  the  molecule  these  two  states  are  coupled  by  an  external 
electric  field,  even  when  spin-orbit  interaction  is  absent.  The  spin-electric  coupling  can  be  viewed 
as  originating  from  a  modified  exchange  constant  5J  induced  by  the  electric  field.  We  find  that 
the  calculated  transition  rate  between  the  chiral  states  yields  an  effective  electric  dipole  moment 
d  =  3.38  x  1CF33C  m  ft  el0_4a,  where  a  is  the  Cu  separation.  For  external  electric  fields  e  ft  108 
V/m  this  value  corresponds  to  a  Rabi  time  r  ft:  1  ns  and  to  a  8J  of  the  order  of  a  few  /reV. 


I.  INTRODUCTION 

Single-molecule  magnets  (SMMs)  have  been  in¬ 
tensively  studied  in  the  last  two  decades  (for  a  re¬ 
view  see  Ref.  Oj) .  At  low  temperature  these  remark¬ 
able  molecules  behave  in  part  like  bulk  magnets 
thanks  to  their  very  long  magnetization  relaxation 
time.  At  the  same  time  SMMs  are  genuine  quan¬ 
tum  systems.  They  display  a  variety  of  non-trivial 
quantum  effects  such  as  the  quantum  tunneling  of 
the  magnetization*^  Berry  phase  interference^  and 
quantum  spin  coherence^.  Due  to  their  double  na¬ 
ture,  SMMs  are  ideal  systems  to  investigate  decoher¬ 
ence  and  the  interplay  between  classical  and  quan¬ 
tum  behavior  £ 

From  the  point  of  view  of  applications,  interest  in 
SMMs  has  been  in  part  spurred  by  the  possibility 
that  these  structures  could  represent  the  ultimate 
molecular-scale  limit  for  magnetic  units  in  high- 
density  magnetic  storage  materials.  More  recently 
SMMs  have  been  recognized  as  promising  building 
blocks  in  molecular  spintronics,  the  emerging  field 
combining  spintronics  and  molecular  electronics^— 
In  particular,  thanks  to  their  long  spin  coherence 
tirne^  SMMs  are  good  candidates  to  realize  spin- 
tronic  devices  that  maintain,  control  and  exploit 
quantum  coherence  of  individual  spin  states.  These 
devices  could  find  important  applications  in  the  field 
of  quantum  information  processing.12,13 

One  key  issue  in  using  SMMs  in  molecular  spin¬ 
tronics  and  quantum  information  processing  is  the 
ability  of  switching  efficiently  between  their  different 
magnetic  states.  The  conventional  way  of  manip¬ 
ulating  magnetic  states  is  by  applying  an  external 
magnetic  field.  However,  this  approach  has  signif¬ 


icant  drawbacks  when  it  comes  to  controlling  mag¬ 
netic  states  at  the  molecular  level.  Quantum  manip¬ 
ulation  of  SMM  requires  application  of  an  external 
field  at  a  very  small  spatial  and  temporal  scale.  It  is, 
however,  very  difficult  to  achieve  such  a  small  scale 
manipulation  using  standard  electron-  spin  control 
techniques  such  as  electron  spin  resonance  (ESR) 
driven  by  ac  magnetic  fields 

One  promising  alternative  to  achieve  control  of 
magnetic  states  at  the  molecular  level  is  to  use  an 
electric  field  instead.  Typically,  by  using  STM  tips 
for  example,  it  is  possible  to  apply  strong  time- 
dependent  electric  fields  in  sub-nano  regions,  with 
time  scales  of  1  nsd4’15  Clearly,  since  electric  fields 
do  not  couple  directly  to  spins,  it  is  essential  to  find 
efficient  mechanisms  for  spin-electric  coupling  as  well 
as  real  SMMs  where  this  mechanism  can  be  at  play. 
In  principle  an  electric  field  can  interact  with  spins 
indirectly  via  the  spin-orbit  interaction.  However, 
since  the  strength  of  the  coupling  scales  like  the  vol¬ 
ume  of  the  system,  this  mechanism  is  not  the  most 
efficient  one  for  manipulating  SMMs. 

Recently,  it  has  been  proposed^  that  in  some 
molecular  antiferromagnets  lacking  inversion  sym¬ 
metry,  such  as  the  triangular  antiferromagnetic 
{Cuz}  and  other  odd-spin  rings,  an  electric  field 
can  efficiently  couple  spin  states  through  a  com¬ 
bination  of  exchange  and  chiralilty  of  the  spin- 
manifold  ground  stated.  The  {C113}  molecule,  while 
largo17’18,  reduces  to  a  simple  model  composed  of 
three  identical  spin  s  =  1/2  Cu  cations  coupled 
by  an  antiferromagnetic  (Heisenberg)  exchange  in¬ 
tegration.  Its  ground  state  consists  of  two  total-spin 
5=1/2  doublets  of  opposite  spin  chirality,  degener¬ 
ate  in  the  absence  of  spin-orbit  interaction.  Accord- 
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ing  to  an  analysis  based  on  group  theory;16!19  due 
to  the  lack  of  inversion  symmetry,  an  electric  field 
can  couple  states  of  opposite  chirality  through  the 
dipole  operator,  even  when  spin-orbit  interaction  is 
absent.  In  the  presence  of  an  additional  small  dc 
magnetic  field  that  mixes  the  spin  states,  this  spin- 
electric  coupling  will  then  generate  efficient  electric 
transitions  from  one  spin  state  to  another. 

An  intuitive  picture  of  this  coupling  is  the  fol¬ 
lowing.  Since  a  spin  S  =  1/2  triangular  antiferro- 
magnet  is  frustrated,  there  exist  three  energetically 
degenerate  antiferromagnetic  spin  configurations  for 
Sz  =  1/2  and  three  for  Sz  =  —1/2.  Both  ground 
state  chiral  eigenstates,  with  a  given  value  of  Sz,  are 
appropriate,  equally  weighted,  linear  combinations 
of  these  three  frustrated  spin  configurations.  Each  of 
these  three  configurations,  if  prepared,  would  have  a 
dipole  moment  with  the  same  magnitude  that  points 
from  the  antiparallcl  sites  to  the  midpoint  between 
the  two  parallel  sites.  While  the  net  dipole  moment 
of  the  two  chiral  eigenstates  is  zero,  the  dipole  tran¬ 
sition  matrix  element  between  them  is  not  and  it  is 
simply  related  to  the  magnitude  of  the  permanent 
dipole  moment  of  the  energetically  degenerate  frus¬ 
trated  configurations. 

In  practice  the  relevance  of  this  spin-electric  mech¬ 
anism  depends  on  the  coupling  strength  of  the  chiral 
states  by  the  electric  field,  i.e.  on  the  value  of  the 
dipole  moment  of  the  frustrated  spin  configurations. 
Theoretically  this  is  an  issue  that  only  a  micro¬ 
scopic  calculation  for  the  specific  molecule  can  ad¬ 
dress.  The  main  objective  of  this  work  is  to  calculate 
the  strength  of  this  coupling  for  the  {C113}  molecule 
using  a&-initio  methods.  Our  approach  is  based 
on  Spin  Density  Functional  Theory  (SDFT),  imple¬ 
mented  in  the  NRLMOL  codes,  which  has  been  very 
successful  in  describing  the  electronics  and  magnetic 
properties  of  Mni2-acetate  and  other  SMMs.— 
Recently  SDFT  implemented  in  NRLMOL  has  been 
used  in  a  first-principle  study  of  quantum  transport 
in  a  Mni2  single-electron  transistor^ 

Our  results  show  that  indeed  the  crucial  electric- 
dipole  moment  is  not  negligible  in  {C113}  and  it 
would  correspond  to  characteristic  Rabi  times  of  1 
ns  in  the  presence  of  typical  electric  fields  generated 
by  STM  tips.  As  originally  suggested  in  Ref.  0. 
the  spin-electric  coupling  can  be  interpreted  as  due 
to  a  modified  exchange  interaction  brought  about 
by  the  electric  field.  Although  here  we  only  address 
the  specific  case  of  {C113},  our  paper  introduces  a 
methodology  that  can  be  followed  in  a  systematic 
study  of  other  SMMs  without  inversion  symmetry. 

The  paper  is  organized  as  follows.  In  section  |TT] 
we  discuss  the  electronic  and  magnetic  properties  of 
triangular  {C1L3}  molecule  based  on  a&-initio  calcu¬ 
lations  and  show  that  the  low-energy  quantum  prop¬ 


erties  of  the  molecule  can  be  described  by  an  effectiv? 
three-spin  s  =  1/2  Heisenberg  model  with  antiferro¬ 
magnetic  coupling.  In  section  III  Bl  we  review  the 
underlying  mechanism  of  spin-electric  coupling  in 
{CU3}  antiferromagnet,  based  on  the  effective  spin 
Hamiltonian.  The  first-principle  computation  of  the 
spin-electric  coupling  and  electric  dipole  moment  of 
{Crta}  is  presented  in  section  IIVI  In  section  IIVBI 
we  discuss  the  effect  of  the  electric  field  on  the  ex¬ 
change  coupling.  Finally  we  present  the  summary  of 
our  work  in  section  El 


II.  ELECTRONIC  AND  MAGNETIC 
PROPERTIES  OF  {Cu3} 

A.  Microscopic  description  of  the  molecule 

The  {CU3}  molecule  that  we  are  interested  in 
has  chemical  composition  Nai2[Cu3(AsWg033)2- 
3H20]-32H2Oii.  This  molecule  has  been  studied 
experimentally  by  different  groups j—  The  three 
Cu2+  cations  form  an  equilateral  triangle  and,  as  we 
show  below,  are  the  sites  of  three  identical  s  =  1/2 
quantum  spins.  The  frontier  electrons  on  each  of 
these  sites  have  primarily  d  character.  The  bridging 
atoms  consist  of  predominantly  paired  electrons  and 
are  only  polarized  to  the  degree  that  the  same-spin 
states  hybridize  with  the  unpaired  d-electrons  on  the 
Cu  sites.  Due  to  the  localized  nature  of  transition- 
metal  3d  states,  direct  exchange  stabilization  due 
to  parallel  neighboring  states  is  expected  to  be  ex¬ 
ponentially  small.  Therefore,  unless  the  frontier  d- 
electrons  are  spatially  orthogonal  by  symmetry  to 
the  d-electrons  on  other  sites,  antiferromagnetic  or¬ 
dering  between  electrons  on  a  pair  of  neighboring 
Cu  atoms  is  energetically  preferred  due  to  the  in¬ 
crease  in  the  system’s  kinetic  energy,  induced  by  or¬ 
thogonality  constraints,  when  neighboring  states  are 
parallel. 

Although  the  spin  model  of  three  exchange- 
coupled  spin  1/2  is  quite  useful  to  understand  the 
magnetic  properties  of  the  {Cus}  SMM,  all  the  other 
atoms  in  the  molecule  are  essential  for  its  geometri¬ 
cal  stability  and  for  the  resulting  superexchange  in¬ 
teraction  among  the  spins  at  the  Cu  sites.  A  proper 
afe-initio  description  of  the  molecule  must  therefore 
include  to  a  certain  extent  all  these  atoms. 

Building  a  suitable  model  of  the  molecule  is  a  con¬ 
siderable  challenge  since  the  model  molecule  should 
preserve  the  essential  physics.  We  have  constructed 
the  molecule  by  preserving  the  D3/j  symmetry  of  the 
polyanionic  part  of  the  molecule  as  observed  in  the 
experiment]17’18  Three  of  the  twelve  Na  atoms  of 
the  molecule  are  placed  at  the  belt  region  of  the 
molecule.  These  three  are  the  most  important  of  all 


the  Na  atoms  for  the  stability  of  the  belt  region  of 
the  molecule.  There  is  some  uncertainty  in  the  posi¬ 
tion  of  the  Na  atoms  but  we  have  placed  eight  of  the 
remaining  nine  Na  atoms  in  a  way  to  preserve  the 
D3/l  symmetry.  The  last  Na  atom  is  replaced  by  a  H 
atom  and  is  placed  at  the  center  of  the  molecule  to 
maintain  the  charge  neutrality  of  the  valance  elec¬ 
trons.  The  model  of  the  molecule  used  in  this  calcu¬ 
lation  is  shown  in  Fig.  [TJ 


FIG.  1:  Model  of  the  {CM3}  molecule  with  chemical 
composition  NanH[Cu3(AsW9  033)2-  3H2O]  used  in  this 
work.  Xcrysden  visualization  tool^  is  used  for  this  fig¬ 
ure. 

We  have  relaxed  the  geometry  using  the  a&-initio 
package  NRLMOL2™i2-  that  uses  a  Gaussian  basis 
set  to  solve  the  Kohn-Sham  equations  using  PBE- 
GGA  approximation^  All-electron  calculations  are 
performed  for  all  elements  of  the  molecule  except 
for  tungsten,  for  which  we  have  used  pseudo  poten¬ 
tials.  The  relaxation  is  first  performed  by  setting  the 
net  total  spin  of  the  molecule  to  S  =  3/2  and  then 
by  changing  the  net  spin  to  1/2.  Self-consistency  is 
reached  when  the  total  energy  is  converged  to  10-6 
Hartree  or  less. 

The  density  of  states  of  the  molecule  is  shown  in 
Fig.  [21  The  HOMO-LUMO  gap  for  the  majority 
spin  is  calculated  to  be  about  0.78  eV  and  that  for 
minority  spin  is  about  0.58  eV.  Although  in  our  cal¬ 
culations  we  have  used  an  equilateral  arrangement  of 
the  three  Cu  atoms,  it  is  found  experimentally  that 
the  {Czt3}  molecule  in  the  ground  state  is  slightly 
distorted  into  an  isosceles  triangle.—  Since  the  cal¬ 
culated  HOMO-LUMO  gap  for  the  equilateral  con¬ 
figuration  is  relatively  large,  the  distortion  is  likely 
to  be  due  to  magnetic  exchange  rather  than  to  the 
Jahn- Teller  effect. 

One  important  result  of  our  calculations,  after  the 


FIG.  2:  Density  of  states  of  {CM3}  molecule.  HOMO- 
LUMO  gaps  for  majority  and  minority  spins  are  shown 
in  the  inset. 


geometry  relaxation  have  been  implemented,  is  that 
the  ground  state  of  the  system  is  anti-ferromagnetic, 
with  a  net  total  spin  of  1/2  in  accordance  with 
experiment.2®  The  ground  state  energy  is  lower  by 
about  5.6  meV  relative  to  spin  S  =  3/2  configura¬ 
tion.  This  allows  us  to  assign  an  exchange  constant 
J  ss  5  meV  to  the  three-site  Heisenberg  spin  model 
mentioned  above  (see  also  next  Section). 

The  calculated  magnetization  density  of  the  re¬ 
laxed  molecule  shows  the  presence  of  three  electron- 
spin  magnetic  moments  ^  ~  0.55/zb,  i  =  1,  2, 3,  es¬ 
sentially  localized  at  the  three  Cu  atom  sites.  Note 
that  the  orbital  moments  are  quenched.  These  re¬ 
sults  confirm  that  the  low-energy  properties  of  the 
{Czz3}  molecule  can  be  approximately  described  by 
an  effective  spin  Hamiltonian  of  three  spins  s  =  1/2 
localized  at  the  Cu  sites. 

The  exchange  coupling  between  two  Cu  atoms  is 
indirect  and  follows  a  superexchange  path^  along 
Cu-O-W-O-W-O-Cu  as  shown  in  Fig.  [3]  -  see  the  ye¬ 
llow  line  connecting  the  atoms.  To  understand  this 
coupling  mechanism  we  focus  on  one  of  the  three 
C11O5  complexes  of  the  molecule  (shown  inside  the 
circle  in  Fig.  0.  Because  of  the  square-pyramidal 
C4v  point-group  symmetry  of  this  complex,  the  dxy, 
dxz,  dyz  states  of  Cu  have  lower  energies  compared 
to  the  dx2_y2  and  dz 2  states.  Moreover,  our  calcu¬ 
lation  shows  that  the  axial  Cu-0  distance  (2.35  A) 
in  each  unit  is  larger  than  the  four  equatorial  Cu- 
O  distances  (1.93  A).  Thus  the  energy  of  dz 2  state 
is  lower  than  dx 2_y2  state  and  the  unpaired  d  elec¬ 
tron  of  the  Cu2+  ion  resides  in  dx2_y2  state  that  is 
directed  along  the  equatorial  Cu-0  vectors.  There¬ 
fore,  the  exchange  coupling  between  two  Cu  atoms 
involves  three  O  atoms  and  two  W  atoms. 

The  magnetic  moment  calculations  of  the  atoms 
of  {Czz3}  molecule  also  support  the  superexchange 
path.  The  magnetic  moments  at  the  O  and  W  atoms 
on  this  path  is  much  smaller  than  at  the  Cu  sites, 


FIG.  3:  (Color  online)  Superexchange  coupling  between 
two  Cu  atoms.  The  yellow  line  connecting  two  Cu  atoms 
through  three  O  and  two  W  atoms  shows  the  path  along 
which  spin  coupling  between  Cu  atoms  is  mediated.  The 
numbers  near  the  atoms  are  the  magnetic  moment  (in 
units  of  /x_b)  of  the  atoms  along  the  exchange  path. 


FIG.  4:  a)  The  three  spin  configurations  of  the  molecule 
associated  with  total  spin  projection  Sz  =  +1/2.  b) 
The  two  chiral  states  formed  from  a)  with  chirality  +1 
and  -1,  respectively. 


but  still  2  order  of  magnitude  larger  than  at  atoms 
not  belonging  to  this  path. 


B.  Effective  spin  Hamiltonian  description 

Based  on  the  results  of  the  a&-initio  calculations, 
the  low-energy  properties  of  the  {CM3}  molecule  can 
be  described  by  the  following  quantum  spin  Hamil¬ 
tonian 

3  3 

H0  =  E  ‘  Sj+l  +  E°  2,2- {-1  *  ®2  X  ,  (1) 
2=1  2=1 

where  J  is  the  exchange  parameter,  D  is  the 
Dzyaloshinski  vector  and  s,  are  three  spins- 1/2, 
located  at  the  Cu  sites.  The  first  term  in  the 
Hamiltonian  is  an  isotropic  Heisenberg  model.  The 
geometry-relaxation  and  electronic-structure  calcu¬ 
lations  showed  that  the  Cu  atoms  form  an  equilat¬ 
eral  triangle  with  a  very  small  intrinsic  deformation. 
Since  the  atomic  environment  around  each  of  the 
the  three  Cu-Cu  bonds  is  the  same,  we  take  the 
three  exchange  constants  J^y+i  to  be  the  same  value, 
J.  On  the  basis  of  the  splitting  between  the  ferro¬ 
magnetic  and  antiferromagnetic  configurations  dis¬ 
cussed  in  the  previous  section,  J  is  positive  and  ~  5 
meV.  The  second  term  in  Eq.  HD  is  the  anisotropic 
Dzyaloshinski-Moriya  exchange  interaction  origina¬ 
ting  from  spin-orbit  interaction.  Its  strength  iD^j+i  | 
is  at  least  one-order  of  magnitude  smaller  than  the 
isotropic  exchange  constant  J,  and  we  will  disregard 
it  for  the  moment. 

The  ground  state  of  Eq.  m  is  total  spin  S'  =  1/2 
manifold,  which  can  be  constructed  in  terms  of  six 
degenerate  spin  configurations,  three  associated  with 


Sz  =  +1/2  and  the  other  three  associated  with 
Sz  =  —1/2.  Fig.  [4]  shows  the  three  possible  spin 
configurations  associated  with  Sz  =  +1/2. 

The  total-spin  S  =  3/2  four-dimensional  subspace 
has  an  energy  of  order  J  above  the  ground-state 
manifold. 

Within  the  S  =  1/2  ground-state  manifold,  we  can 
construct  two  degenerate,  linearly  independent  dou¬ 
blets.  Specifically  the  two  Sz  =  +1/2  states  (shown 
in  Fig.  IHb))  are 

(£’+,+-)  =  +w2’J'ti+]  , 

|-E-,+-)  =  -j=  [’Fj.ft  +  +  wtH-i-tf.]  ,  (2) 

while  the  Sz  =  —1/2  states  are 

I-E+j-j)  =  +W^4.U  +  > 

=  ^=[^f44.  +  w2'F1/|q  +  WF144-]  ,  (3) 

where  w  =  e1?.  The  quantum  numbers  E+  and 
E_  specify  the  so  called  handness  or  chirality  of  the 
states  |  E±,  M ),  which  are  eigenstates  of  the  chirality 
operator 


Cz  =  -^=si  •  s2  x  s3  , 


(4) 


with  eigenvalues  ±1  respectively.  It  is  useful  to  in¬ 
troduce  also  the  other  two  components  of  the  chiral 
vector  operator 

2 

Cx  =  —  g (si  •  s2  -  2s2  •  s3  +  s3  ■  si)  ,  (5) 


Cy  = -^(Si  •  s2  -  s3  •  si)  ,  (6) 


and  the  ladder  operators  C±  =  Cx  ±  iCy.  Note 
that  [Ci,Cm]  =  i2eimuCn  and  [ChSm\  =  0.  Here 
£imn  is  the  Levi-Civita  symbol.  The  ladder  opera¬ 
tors  reverse  the  chirality  of  the  states:  C±\Ezp,  M)  = 
\E±,M).  Thus  C  behaves  exactly  like  the  operator 
S  (for  S  =  1/2)  in  chiral  space. 

In  the  microscopic  description  of  the  molecule  im¬ 
plemented  within  density  functional  theory  via  the 
NRLMOL  code,  the  chiral  states  defined  in  Eqs.  © 
and  have  to  be  understood  as  being  composed 
both  of  a  spin  and  an  orbital  part. 

We  conclude  this  section  with  an  observation  of 
the  DM  interaction.  As  shown  in  Ref.  Util.  the  DM  in¬ 
teraction  within  the  5  =  1/2  ground  state  manifold 
takes  the  simple  form  £?dm  =  Ago  CZSZ,  where  Ago 
is  the  effective  spin  orbit  coupling  constant.  Thus 
equal-spin  states  of  opposite  chirality  are  split  by 
2Ago- 


III.  SPIN-ELECTRIC  EFFECT  IN  {Cu3} 

A.  Absence  of  inversion  symmetry  and 
coupling  of  ground-state  chiral  states 


turns  out  that  in  triangular  spin-1/2  antiferroma^ 
nets  the  coupling  is  non-zero.  On  the  other  hand 
in  pentagon  spin  1/2  antiferromagnets,  the  coupling 
vanishes,  unless  spin-orbit  interaction  is  included^ 
We  focus  now  on  the  spin-electric  coupling  of  chi¬ 
ral  states  in  {Czi3}.  In  the  presence  of  an  external 
electric  field  e,  The  Hamiltonian  acquires  the  addi¬ 
tional  electric-dipole  term  He  =  )T/:  er j  •  e  =  eR  •  e. 
where  e  is  the  electron  charge  and  r j  is  the  coordi¬ 
nate  of  the  ith  electron. 

In  the  subspace  of  spin  projection  Sz  =  1/2  of  the 
ground-state  manifold,  which  is  invariant  for  He  the 
perturbed  Hamiltonian  Hq  +  He  can  be  expressed  in 
the  basis  of  the  chiral  states  as 

H  =  H o  ±  Hs 

{E+,  +\\Hq\E+,  +i)  (E+,  +^\He\E_,  ±|) 
(U_,+i|Re|S+,+i)  (.E_,-4|flo|£L,+i) 

(7) 

A  similar  expression  holds  for  the  Sz  =  —1/2  sub¬ 
space.  The  eigenvalues  of  H  are 

££(e)  =  i£(0)±|d-e|,  (8) 

2  2 

with  Ef( 0)  =  (E±,+^\Hq\E±,+^),  and  the  corre¬ 
sponding  eigenstates 


The  triangular  spin-1/2  antiferromagnet  {Czi3} 
belongs  to  the  class  of  antiferromagnetic  rings  with 
an  odd  number  of  half-integer  spins.  In  these  sys¬ 
tems,  the  lack  of  inversion  symmetry  of  the  molecule 
as  a  whole  implies  that  the  ground-state  is  a  four¬ 
dimensional  manifold,  whose  basis  states  | E±,SZ  = 
±1/2)  are  characterized  by  the  spin  projection  Sz  = 
±1/2  and  by  the  chirality  Cz  =  ±1  (which  we  also 
label  as  E±  ).  In  contrast,  antiferromagnetic  rings 
with  an  even  number  of  spins  have  non-degenerate 
5  =  0  singlet  ground  state.  According  to  the  origi¬ 
nal  proposal  in  Ref.  [Ml  in  odd-spin  rings  the  two 
states  of  opposite  chirality  \E±,SZ  =  M)  can  be 
coupled  linearly  by  an  external  electric  held,  even 
in  the  absence  of  spin-orbit  interaction.  In  order 
for  electric  coupling  to  be  non-zero,  other  criteria 
must  be  satisfied*^  First  of  all,  permanent  electric 
dipoles  dij  must  be  present  on  the  bridges  that  me¬ 
diate  the  coupling  of  spin  s,  and  s j.  A  necessary 
(although  not  sufficient)  condition  for  this  is  that 
the  superexchange  bridge  that  magnetically  couples 
Si  and  Sj  lacks  a  center  of  inversion  symmetry.  Even 
when  local  dipole  moments  are  present  on  individ¬ 
ual  bridges,  the  resulting  final  spin-electric  coupling 
between  chiral  states  depends  in  a  nontrivial  way 
on  the  overall  symmetry  of  the  molecule.  The  best 
way  to  settle  this  issue  is  to  carry  out  a  system¬ 
atic  symmetry  analysis  based  on  group  theory.  It 


xf(e)j  = 


V2 


'^■+5>±JST^-+5> 


(9) 

Here  we  have  introduced  the  electric  dipole  matrix 
element  d,  which  couples  states  of  opposite  chirality 
(but  with  the  same  spin  projection) 


d  —  (E+,  ±-|eR|I?-,  +— )  .  (10) 


For  the  specific  example  of  {C1L3}  molecule  only 
the  matrix  elements  of  X  and  Y  components  of  R 
are  nonzero  and 

(-®+>  +^leA!|-E-,  ±^)  =  i{E+,  +-\eY\E_,  +-)  =  -j=  , 

(11) 

where  d,  =  |d|. 

The  matrix  element  in  Eq.  m  is  the  key  quantity 
in  the  spin-electric  coupling  mechanism.  Substitut¬ 
ing  the  expressions  for  the  chiral  states  from  Eqs.  © 
and  using  the  orthogonality  of  spin  states  we  obtain 

d  =  g  ( ( ^  Ttt  |  eR.  |  4_tt)  +  + 

w2(^m|eR|^m))  .  (12) 


Evaluating  the  dipole  matrix  element  between  two 
states  of  opposite  chirality  is  therefore  equivalent  to 


calculating  the  dipole  moment  of  each  of  the  three 
spin  configurations.  This  matrix  element  determines 
the  strength  of  spin-electric  coupling  and  we  are  pri¬ 
marily  interested  in  calculating  this  quantity  by  ab- 
initio  methods. 

Finally,  note  that  all  the  matrix  elements  of 
the  electric  dipole  operator  eR  are  identically 
zero  in  the  S  =  3/2  subspace.  This  is  obvi¬ 
ous  since  (’F^-|-|eR|'I'^y[-)  and  3 ((^4.ttleR|^,4.tt)  + 
('Ft4.t|eR|’Fnt}  +  (^m|eR|Tm))  are  both  zero  by 
symmetry.  We  will  confirm  this  result  by  direct  ab- 
initio  calculations. 


B.  Effective  spin  Hamiltonian  description 

The  effect  of  the  electric  Held  on  the  the  low- 
energy  spectrum  of  {CM3}  can  be  recast  in  the  form 
of  the  effective  spin  model  introduced  in  Sec.  mu 
Since  the  electric  dipole  operator  has  nonzero  matrix 
elements  only  in  the  ground-state  manifold,  where 
it  couples  states  with  equal  spin  components  and 
opposite  chirality,  we  expect  that  the  spin-electric 
Hamiltonian  He  can  be  rewritten  as  a  linear  com¬ 
bination  of  the  ladder  operators  C± .  By  comparing 
the  matrix  elements  of  Hs  given  in  Eq.  (ITT)])  and  (fill) 
with  the  action  of  C±  on  the  chiral  states,  one  can 
show  thatiS 

=  11  ’  (13) 

where  e'  =  Rz((/))(7tt/6  —  2 9)e,  with  R(<j>)  being  the 
matrix  representing  a  rotation  by  an  angle  <fi  around 
the  2-axis,  and  0  being  the  angle  between  the  in¬ 
plane  component  £y  of  the  electric  field  and  the  bond 
Si  —  S2.  By  using  Eq.  (O  and  J6|  we  can  now  rewrite 
C||  =  ( CXlCy )  in  term  of  spin-operators  s,;  and  we 
obtain— 

3 

Hls  =  ^2,5Jii+1{e)si  •  s j  ,  (14) 

i 

where  the  modified  exchange  parameters  take  the 
formic 

4-/7  27r 

8Ja+i(e)  =  ^y=|£|||cos(— i  +  9)  .  (15) 

This  expression  of  the  effective  electric-dipole 
Hamiltonian  suggests  a  transparent  physical  in¬ 
terpretation  of  the  spin-electric  couping  mecha¬ 
nism.—  An  external  electric  field  changes  the 
charge  distribution  of  the  {CM3}  molecule  which, 
in  turn,  changes  the  exchange  interaction  between 
neighboring  atoms.  Since  the  modified  exchange  in¬ 
teraction  does  not  commute  with  Hq,  it  can  cause 


transitions  between  chiral  states  within  the  ground 
state  manifold. 

In  Eq.  (HU),  £n  is  the  projection  of  electric  field 
on  the  CM3  plane  (in  our  case  ey  =  e),  i=l  and 
6  =  30°  is  the  angle  between  e  and  the  line  joining 
Cmi  and  Cu2 ■  Finally,  note  that  Eqs.  m  and  m 
provide  an  estimate  of  the  dependence  of  the  ground- 
state  energy  as  function  of  the  electric  field.  Since 
in  the  absence  of  spin-orbit  coupling  the  electric- 
dipole  Hamiltonian  has  zero  matrix  elements  in  the 
S  =  3/2  subspace,  Eq.  (fl5l)  gives  us  an  estimate  of 
the  dependence  of  the  exchange  constant  J  (propor¬ 
tional  to  the  splitting  between  the  S'  =  1/2  ground 
state  and  S  =  3/2  excited  state)  on  e. 


IV.  AB-INITIO  EVALUATION  OF  THE 
SPIN-ELECTRIC  COUPLING 

A.  Calculation  of  the  electric  dipole  moment 

To  construct  the  chiral  states  of  the  full  {CM3} 
molecule,  we  have  calculated  the  ground  state  of  the 
molecule  for  different  spin  configurations,  as  shown 
in  Fig.  [I]  Although  there  are  two  doublets  of  chiral 
states  for  the  triangular  arrangement  of  three  spin 
1/2  atoms,  in  this  calculation  we  have  used  only  one 
doublet  associated  with  the  spin  projection  +1/2, 
since  we  are  interested  in  coupling  between  states  of 
opposite  chirality  with  the  same  spin  projection. 

To  study  the  spin-electric  effect  we  have  applied 
an  external  field  along  the  perpendicular  bisectors 
between  positions  2  and  3  of  the  CM3  triangle  shown 
in  Fig.  0  and  have  calculated  the  corresponding 
ground  state  energy  self-consistently  for  different 
spin  configurations.  We  have  kept  the  direction  of 
the  field  relative  to  coordinate  axes  fixed,  and  have 
changed  the  orientation  of  the  spins  at  the  Cu  atoms 
to  generate  the  three  possible  spin  configurations  of 
the  {CM3}  molecule. 


FIG.  5:  The  direction  of  the  applied  electric  fields  used 
in  this  calculation. 


Our  calculations  show  that  {CM3}  molecule  in 
the  spin  Sz  =  3/2  state  does  not  have  any  per¬ 
manent  electric-dipole  moment.  On  the  other  hand 
each  of  the  three  frustrated  spin  Sz  =  1/2  con¬ 
figurations  have  a  small  permanent  (i.e,  zero-field) 
dipole  moment,  as  expected  from  the  general  dis¬ 
cussion  of  Sec.  III.  The  three  moments  have  all  the 
same  magnitude  but  their  directions  are  along  the 
perpendicular  bisector  of  C113  triangle  and  between 
two  Cu  atoms  with  parallel  spin  alignments.  The 
relative  orientations  of  these  moments  along  with 
components  are  shown  in  Fig.  [G]  The  fact  that  the 
Sz  =  3/2  state  does  not  have  permanent  dipole  mo¬ 
ment  whereas  Sz  =  1/2  states  do,  suggests  that  the 
dipole  moments  are  solely  due  to  spin  effects. 


FIG.  6:  Dipole  moments  of  three  spin  configurations 
and  their  relative  angles.  Pprt  =  ( ^ I  eFt |  j.-p-r ) , 

Put  =  (®utleRl®t4.t)  and  Pm  =  (^mkR-l^m)  are 
the  moments  corresponding  to  the  spin  configurations  of 
Fig.  |4}i- 

In  the  presence  of  an  electric  field  the  energies 
of  the  {CM3}  molecule  are  slightly  lower  when  field 
is  between  two  Cu  atoms  with  parallel  spins  than 
for  the  other  two  spin  configurations,  where  the 
field  is  between  two  Cu  atoms  with  anti-parallel 
spin  alignments.  This  difference  in  energy  is  due  to 
the  direction  of  permanent  moment  relative  to  the 
induced  moment.  We  have  calculated  the  perma¬ 
nent  dipole-moment  of  the  ground  state  spin  config¬ 
uration  by  fitting  the  dependence  of  energy  of  one 
of  the  Sz  =  1/2  spin  configurations  with  external 
field,  as  shown  in  Fig.  [7]  The  calculated  values 
of  the  permanent  dipole  moment  and  polarizability 
of  {Czt3}  molecule  are  p  =  4.77  x  10-33C.m  and 
a  =  1.025  x  10 ~38C.m2 /V ,  respectively.  Although 
there  is  no  experimental  value  of  polarizability  avail¬ 
able  for  {CM3},  polarizabilities  within  DFT  calcula¬ 
tions  are  generally  accurate  to  1-3  percent. 

The  value  p  extracted  from  this  fitting  is  consis¬ 


tent  with  the  direct  calculation  of  the  electric  dipol? 
moment  of  the  three  spin  configurations  at  zero  field, 
implemented  in  the  NRLMOL. 
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FIG.  7:  Electric  field  dependence  of  the  energy  for  one 
of  the  three  spin  Sz  =  1/2  spin  configurations.  The  plot 
for  the  other  two  configurations  is  very  similar  and  the 
fitting  yields  essentially  the  same  values  of  p  and  a. 

To  calculate  the  matrix  element  d  given  in 
Eq.  (ED,  we  substitute  the  components  of  the  mo¬ 
ments  for  the  different  spin  configurations  of  Fig.  [Gj 

d  =  -p[(l  +  oj  cos  p  +  cm2  cos  2/3)x 

O 

+(w  sin  f3  +  to2  sin  2/3)y] 

=  |(x  +  *y).  (16) 

The  magnitude  of  the  dipole  coupling  in  {CM3} 
molecule  is,  therefore 

d  =  -y=  =  3.38  x  10"33C.m  .  (17) 

v  2 

The  efficiency  of  the  {CM3}  molecule  as  a  switch¬ 
ing  device  depends  on  how  fast  an  electric  field  can 
generate  transitions  from  one  chiral  state  to  the 
other.  The  characteristic  (Rabi)  time  for  transitions 
between  the  two  chiral  states  is  given  by 


Here,  h  is  Planck  constant,  d  is  the  dipole  matrix 
element  between  states  of  different  chirality  given  by 
Eq.  (flTl).  and  e  is  the  external  electric  field.  Fig.  [8] 
shows  the  dependence  of  the  Rabi  time  on  external 
field,  with  the  maximum  value  of  ~  50  ns  for  a  field 
£  ~  5  x  106  V/m.  For  larger  fields  of  the  order 
of  ~  108  V/m,  easily  attainable  in  the  vicinity  of  a 
STM  tip,  the  Rabi  time  is  of  the  order  of  1  ns,  which 
is  considered  to  be  a  relatively  fast  control-time  in 
quantum  information  processing. 


FIG.  8:  Electric  field  dependence  of  the  Rabi  time  for 
quantum  transitions  between  the  two  (ground-state. 


B.  Modification  of  the  exchange  coupling  in  an 
electric  field 

To  calculate  the  dependence  of  the  exchange  cou¬ 
pling  J  on  the  electric  field,  we  need  to  determine 
how  the  spin  5  =  1/2  ground  state  and  the  spin 
5  =  3/2  excited  state  depend  on  the  field.  We  de¬ 
fine  the  exchange  energy  J(e)  as  the  difference 

J (e)  =  Ei  (e)  —  ET  (e)  ,  (19) 
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where  E  i  (e)  and  E j  (e)  are  the  energies  of  the  5  = 
1/2  ground-state  and  of  the  spin  5  =  3/2  excited 
state  respectively  in  the  presence  of  an  electric  field. 

Based  on  on  our  discussion  of  Sec.  EUA  [see 
Eq.  ©],  the  energy  of  the  5  =  1/2  chiral  ground- 
state  manifold  and  the  5  =  3/2  excited  state  vs  £ 
are  shown  schematically  in  Fig.  [9]  where  we  have  dis¬ 
regarded  the  quadratic  dependence  of  both  E j  and 
Ei  on  the  field  due  to  the  induced  electric  dipole 
moment. 


FIG.  9:  Schematic  electric-field  dependence  of  the  ener¬ 
gies  of  the  S  =  1/2  chiral  states  and  spin  5  =  3/2  excited 
state,  and  the  exchange  energy  J  defined  in  Eq.  m- 

The  calculation  of  the  electric-field-modified  ex¬ 
change  parameter  using  first-principle  methods  is 


not  completely  straightforward,  since  the  SDFT  ca@ 
dilations  done  within  NRLMOL  allow  us  to  calcu¬ 
late  the  energy  of  a  given  spin  configuration,  whereas 
the  (chiral)  ground-state  is  a  linear  combination  of 
three  possible  spin  configurations.  However,  we  can 
get  an  estimate  of  the  dependence  of  J  on  e  by  ap¬ 
proximating 

Ei  (e)  ~  a?(^4.tt|i5DFT(e)|^,j.tt) 

2 

+a2(^'t4-t \Hdft{£) 

+al(^tt4. \Hdft(e) 

The  coefficients  a's  can  be  obtained  by  expanding 
Xl(e))  in  Eq.  j9])  in  terms  of  the  spin  configura¬ 
tions,  which  leads  to 

Xl(e))  =  +  (w  — 

+(w2  -  wr)%;] 

=  +  Q*2^t4.t  +  ,  (20) 


where  =  -E(l  —  j),  for  the  given  choice  of 

the  electric  field  direction. 

Therefore, 


The  energies  E^  and  Ef\ ^  are  the  same  because 
of  symmetry.  Since  the  difference  between  E^  and 
E-^tf  is  very  small  and  near  the  accuracy  limit  of  our 
calculations,  we  further  approximate  E ~  A'tJ.t- 
The  exchange  parameter  J  becomes 

J(e)  w  Em(e)  -  Em(£)  ,  (22) 

with  E|-|-|-(e)  =  E|((e). 

In  Fig.  [TO]  we  plot  the  electric-field-induced  vari¬ 
ation  of  the  exchange  energy  5J(e)  =  J(e)  —  J(0) 
vs.  e.  The  result  for  SJ  obtained  by  evaluating 
Eq.  (|22D  with  SDFT  is  shown  by  the  red  curve.  For 
this  part  of  the  calculations  the  convergence  crite¬ 
rion  has  been  increased  up  to  10-8  Hartree.  We 
can  see  that  the  dependence  of  J  on  electric  field 
is  quite  small,  and  5J  is  in  the  p,eV  range  for  elec¬ 
tric  fields  £  =  (1  —  10)  x  10'  V/m.  These  energies 
are  not  far  form  the  accuracy  limit  of  our  numerical 
calculations,  which  is  the  reason  of  the  fluctuations 
seen  in  the  plot.  Nevertheless  the  overall  trend  is 


an  increase  of  SJ(s)  with  e,  which  is  approximately 
linear  at  low  fields.  Note  that  the  SDFT  evaluations 
of  S-fttfe)  and  E. contain  a  quadratic  contribu¬ 
tion  in  e  but  this  nearly  cancels  at  small  fields  when 
computing  SJ,  and  it  becomes  appreciable  only  at 
e  >  5  x  107  V/m. 

The  blue  line  in  Fig.  [10]  shows  the  dependence  of 
SJ  on  £  given  by  the  prefactor  of  the  cosine  func¬ 
tion  in  Eq.  m,  which  was  derived  within  the  spin 
Hamiltonian  formalism.  When  plotting  Eq.  m  we 
have  used  the  value  of  d  extracted  from  our  first- 
principle  calculations.  Comparing  the  two  curves, 
we  note  that,  apart  form  the  fluctuations  in  the  nu¬ 
merical  result  mentioned  above,  the  theoretical  and 
numerical  values  for  SJ  are  consistent,  and  both  pro¬ 
cedures  predict  an  overall  increase  of  6  J  with  electric 
field. 


FIG.  10:  Electric-field  dependence  of  the  variation  of 
the  exchange  energy  5J{e)  =  J(e)  —  J(0)  induced  by  the 
field.  The  red  curve  is  the  first-principle  result  obtained 
by  evaluating  Eq.  ([22])  and  the  dashed  black  curve  is 
the  quadratic  fit  of  5J{e).  The  blue  curve  is  a  plot  of 
Eq.  m  with  the  numerical  value  of  d  extracted  from 
the  first-principle  calculations. 


V.  SUMMARY 

In  this  paper  we  have  carried  out  a  first- 
principle  study  of  the  spin-electric  coupling  in  single¬ 
molecule  magnets  (SMMs)  without  inversion  sym¬ 
metry.  Specifically,  we  have  analyzed  the  clear-cut 
case  of  the  {C113}  triangular  antiferromagnet  where, 
because  of  spin  frustration,  the  ground-state  con¬ 
sists  of  two  generate  spin  1/2  doublets  of  opposite 
chirality.  Theory  predicts16:19  that  an  electric  field 
can  couple  these  states,  even  when  spin-orbit  inter¬ 
action  is  absent.  The  main  goal  of  our  work  has  been 
to  compute  how  strong  this  coupling  is. 

Our  calculations  of  the  electronic  structure  of  the 
{C113}  molecule  show  that  the  spin  magnetic  mo¬ 


ments  are  localized  at  the  three  Cu  atom  sites  of  the 
molecule.  The  magnetic  properties  of  the  molecule 
are  correctly  described  by  a  triangular  spin  s  =  1/2 
Heisenberg  antiferromagnet,  with  an  exchange  cou¬ 
pling  J  of  the  order  of  5  meV  that  separates  the 
energies  of  the  spin- S'  =  1/2  ground-state  many- 
fold  and  the  spin-S  =  3/2  excited  states.  In  agree¬ 
ment  with  theoretical  predictions)16’19  we  find  that 
an  electric  field  couples  the  two  ground-state  dou¬ 
blets  of  opposite  chirality,  even  when  spin-orbit  in¬ 
teraction  is  absent.  The  strength  of  the  coupling 
is  linear  in  the  field  and  proportional  to  the  perma¬ 
nent  electric  dipole  moment  d  of  the  three  frustrated 
spin  configurations.  The  calculations  yield  a  value 
of  d  ~  4  x  10”33C  m  ft:  el0~4a  for  {CU3},  where  a  is 
the  Cu  atom  separation.  Corresponding  Rabi  times 
for  electric-field-induced  transitions  between  chiral 
states  can  be  as  short  as  1  ns,  for  electric  fields  of 
the  order  of  10s  V/m,  which  are  easily  produced 
by  a  nearby  STM  tip.  Thus  this  spin-electric  cou¬ 
pling  mechanism  is  of  potential  interest  for  the  use 
of  single-molecule  magnets  in  quantum  information 
processing  as  fast  switching  devices. 

Our  calculations  also  indicate  that  the  presence 
of  an  external  electric  field  modifies  the  exchange 
constant  J.  Typically  the  electric  field  increases  J, 
although  the  energy  scale  of  this  change  is  in  the 
/LteV  range  for  typical  STM-generated  electric  fields. 
Thus  for  this  specific  antiferromagnetic  SMM,  the 
electric  field  cannot  trigger  directly  a  level  crossing 
between  magnetic  states  with  different  total  spin,  as 
suggested  recently  for  other  SMMs;30’31 

This  work  shows  that  a  microscopic  investiga¬ 
tion  of  the  spin-electric  coupling  using  the  NRL- 
MOL  first-principle  code  is  feasible,  and  can  sys¬ 
tematically  implemented  for  a  large  class  of  SMMs 
which  lack  inversion  symmetry.  In  this  paper  we 
have  disregarded  the  effect  of  spin-orbit  interaction 
and  external  magnetic  field.  The  spin-orbit  inter¬ 
action  strength  is  small  compared  to  the  exchange 
coupling  J.  In  the  case  of  {Czi3}  it  simply  intro¬ 
duces  a  small  splitting  between  the  chiral  states,  but 
is  not  expected  to  influence  significantly  the  spin- 
electric  coupling.  However  in  other  antiferromag¬ 
netic  rings  with  an  odd  number  of  spins  spin-orbit 
interaction  is  essential  for  the  very  existence  of  the 
coupling  mechanism^.  Work  to  include  both  spin- 
orbit  interactions  and  an  external  magnetic  field  is 
in  progress.  Together  with  the  group-theory  analysis 
presented  in  Ref.  these  studies  will  be  a  consider¬ 
able  help  in  guiding  future  experiments  and  selecting 
the  most  promising  SMMs  for  applications  in  quan¬ 
tum  information  processing  and  nanospintronics. 
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